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Abstract 

In this paper we explore the constraints from B-physics observables in SUSY mod- 
els of Minimal Flavor Violation, in the large tan (3 regime, for both low and high 
scale supersymmetry breaking scenarios. We find that the rare B-decays b —* sj and 
B s — > can be quite sensitive to the scale M at which supersymmetry breaking is 

communicated to the visible sector. In the case of high scale supersymmetry breaking, 
we show that the additional gluino contribution to the b — > S7 and B s — > /i + fj,~ rare 
decay rates can be significant for large tan/3, \i and M3. The constraints on B u — > tv 
are relatively insensitive to the precise scale of M. We also consider the additional 
constraints from the present direct Higgs searches at the Tevatron in the inclusive 
H/A — * tt channel, and the latest CDMS direct dark matter detection experiments. 
We find that altogether the constraints from B-physics, Higgs physics and direct dark 
matter searches can be extremely powerful in probing regions of SUSY parameter space 
for low Ma and large tan (5, leading to a preference for models with a lightest CP-even 
Higgs mass close to the current experimental limit. We find interesting regions of pa- 
rameter space that satisfy all constraints and can be probed by Higgs searches at the 
Tevatron and the LHC and by direct dark matter searches in the near future. 



1 Introduction 



The next few years promise to be extremely exciting for High Energy Physics because of 
new results coming from the Tevatron collider, the expected start of the LHC and a number 
of dark matter detection experiments. It is hoped that all this experimental data will shed 
some light on the mechanism of electroweak symmetry breaking and possibly on the origin 
of dark matter in the universe. 

Theoretically, one of the more promising scenarios that can explain both questions is that 
of low energy supersymmetry. In particular, the minimal supersymmetric extension of the 
Standard Model (MSSM) with R-parity can both stabilize the electroweak scale and provide 
a cold dark matter candidate (i.e. the lightest neutralino) with a relic abundance that is in 
good agreement with the WMAP value [1] 

n C D M h 2 = 0.105i|];{™ (1) 

However, like most extensions of the Standard Model, the MSSM is highly constrained by 
flavor changing effects, in particular through B-physics observables. These constraints can be 
naturally satisfied if the SUSY breaking terms are approximately flavor diagonal at the scale 
M, at which supersymmetry breaking is communicated to the visible sector, and all flavor 
changing effects are loop induced and proportional to the elements of the CKM matrix of 
the Standard Model. Such supersymmetric extensions of the Standard Model are generically 
called Minimal Flavor Violating (MFV) and have been extensively studied in Refs. [2]-[16j. 
In particular, in Ref. [TBJ the impact of maximal CP-violation and minimal flavor violating 
MSSM is considered. At large tan /3, the ratio of the two Higgs vacuum expectation values 
in the MSSM, flavor changing neutral currents (FCNCs) are induced by the Higgs sector 
through loop effects that can lead to significant deviations in B-physics observables from 
their Standard Model predictions. The B-factories, Belle and Babar, and the Tevatron have 
measured many of these observables and these data put strong constraints on the allowed 
MSSM parameter space. 

Simultaneously, experiments are also trying to discover the footprint of supersymmetry 
through dark matter searches of a stable neutralino. These searches for dark matter have 
also begun to put significant constraints on supersymmetric models by providing limits on 
the spin-independent scattering cross-sections of the lightest neutralino with nuclei. In the 
MSSM, the couplings of the down type quarks to the non-standard Higgs bosons are tan (3 
enhanced. Therefore the t-channel Higgs boson contribution to this cross-section can be 
sufficiently enhanced for small enough values of the non-standard CP-even Higgs boson 
mass and large tan/3. In addition the spin-independent cross-section also depends on the 
size of the Higgsino component of the lightest neutralino which is governed by the Higgsino 
mass parameter fi. The impact of direct dark matter searches on Higgs physics has been 
analysed in Refs. PIHCEHICEH]. 

In this article we study the effect of varying the scale M, at which supersymmetry 
breaking is communicated to the visible sector, on B-physics observables in the context of 
Minimal Flavor Violation. We concentrate on two scenarios of supersymmetry breaking: 
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low-scale, M ~ Msusy, and high-scale, M ~ Mqut, SUSY breaking, where Msusy and 
Mqut represent the scale of the supersymmetric particle masses and the Grand Unification 
scale, respectively. In the case of low scale supersymmetry breaking the flavor changing 
effects are governed by loop induced Higgs mediated currents. In the case of high scale 
SUSY breaking, the soft squark mass terms are logarithmically sensitive to the scale M, 
due to their RG evolution. In particular for large supersymmetry breaking scales the soft 
squark mass parameters pick up off-diagonal contributions proportional to the CKM matrix 
elements. Hence the squark and quark mass matrices cannot be diagonalized simultaneously. 
This mismatch between the quark and squark mass bases induces flavor violating quark- 
squark-gluino couplings that are proportional to the CKM matrix elements, which lead to 
important gluino contributions to both the B s — > yU + and b — > 37 rare decays, in addition 
to those already present when M ~ Msusy- 

In order to analyze the size of the possible gluino effects, we shall study scenarios that 
parametrize the possible flavor violation effects in models of Minimal Flavor Violation with 
a small messenger scale M, of the order of the weak scale, and with a large scale M, of the 
order of the GUT scale, respectively. In the first scenario, we shall assume no flavor violating 
quark-squark-gluino couplings. In the second scenario, we shall assume a left-handed squark 
mass matrix that is diagonalized together with the up Yukawa coupling matrix, as would be 
the case if the down Yukawa effects in the RG evolution of the soft masses were neglected 
compared to those of the up Yukawas. The effect of the non-diagonal left-handed down 
type quark-squark-gluino vertices on the B s — > fi + fi~ rare decay within this approximation 
has been previously computed in Ref. [7j. In this article we derive an analytic formula for 
the gluino contribution to the b — > 57 rare decay for large values of tan (3, within the same 
approximation. The validity of this approximation will be discussed in section 12.1.21 

In addition we also study the interplay between the B-physics constraints from the B u — > 
TV, B s — > and the b — > 57 rare decays and the recent direct dark matter detection limits 

from CDMS [20]. Let us stress here that low energy SUSY breaking scenarios lead to a light 
gravitino and therefore the CDMS constraints would not apply. We find that combining the 
limits from B-physics observables, dark matter detection experiments at CDMS and inclusive 
H/A — > tt searches at the Tevatron [24] yields interesting constraints on the Ma — tan/? 
and X t — fi plane, where M^ is the CP-odd Higgs mass and X t is the stop left-right mixing 
parameter. We find regions of parameter space that satisfy all these constraints and can be 
probed by Higgs searches at the Tevatron and by direct dark matter searches in the near 
future. 

The paper is organized as follows, in Section [2] we consider the effect of the scale M 
on B-physics observables in minimal flavor violating MSSM. In particular we present the 
additional gluino contributions to the B s — > fi + fi~ and b — > 57 rare decays. The complete 
calculation of the gluino contribution to the b —>■ sj decay can be found in Appendix [A] 
We also give a brief theoretical overview of the relevant direct dark matter detection cross- 
section. In Section [3] we consider different parametric scenarios that can satisfy the B- 
physics experimental constraints, the limits coming from inclusive H/A — > tt searches at 
the Tevatron and the direct dark matter detection limits from CDMS. In particular we 
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explore the dependence of our results on the scale at which supersymmetry breaking is 
communicated to the visible sector. In Section H] we present our conclusions. 



2 Basic Theoretical Setup 

2.1 B-physics Constraints and messenger mass scale M 

The FCNCs induced by loops of squarks depend on the flavor structure of the soft squark 
mass parameters which, in MFV, is closely tied to the scale at which supersymmetry breaking 
is communicated to the visible sector 0. Assuming the squark masses are flavor independent 
at high scales, the only one-loop corrections that violate flavor are due to RG effects governed 
by the up and down Yukawa matrices, since the gauge interactions are flavor blind. The 
corrections to the left-handed soft SUSY breaking mass parameter at one-loop are given 
by 



AMi 



8tt 



1 ' Mj + M? R + M 2 Hu (0) + A 2 ) Y^Y U + 
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Ml + Ml + Ml( ) + Al) YjY d ] log (-^-) , (2) 

where M? denotes the left-handed squark mass matrix, M? R (M"~ ) is the right-handed up 
(down) squark mass matrix, Mfj (0) and A are the Higgs soft supersymmetry breaking 
and squark-Higgs trilinear mass parameters, respectively, at the messenger scale M, at which 
supersymmetry breaking is transmitted to the observable sector, and Msusy is the charac- 
teristic low energy squark mass scale. Similarly, the right-handed up and down squark mass 
matrices, receive one-loop Yukawa- induced corrections proportional to 

ML = -5? ( M I + <• + M "- (0) + A >) y » y - log (m^) • (3 > 

and 

respectively. Hence the corrections to the right-handed soft mass parameters are diagonal 
in the quark basis, but the left-handed soft mass parameters of the down squarks pick up 
off-diagonal contributions proportional to the CKM matrix elements. The size of these 
corrections depends on the scale M at which SUSY breaking is communicated to the visible 
sector. If M is of the order of Msusy then these corrections are small and if M ~ Mqut 
then these corrections can be substantial. In this section we consider the effect of these two 
scenarios on three B-physics processes b — > sj, B u — > tv and B s — > n + ji~ . 



Unlike Ref. , we are considering the case where effects from the hidden sector are small 
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2.1.1 M ~ Msusy 

In the case M ~ Msusy the squark mass matrices are approximately block diagonal which 
leads to all the neutral Higgs induced FCNCs being proportional to the chargino-stop loop 
factor h 2 ey, with [5] 

ey ~ Yg^^Co(m|,m|, \fi\ 2 ) (5) 

where 

C Q (x,y,z) = rlog(y/x) + - -\og{z/x). (6) 

(x-y)(z-y) (x-z)(y-z) 

No flavor changing effects are produced by contributions from the gluino down squark loop 
as they are purely flavor diagonal 

4 - ^M3^o(m^,mJ /2 ,|M 3 | 2 ), (7) 

where i and mj 2 are the 1 th down squark mass eigenstates. The effective flavor changing 
strange-bottom-neutral-Higgs coupling is [5j [8j [16] 



,s .:>,2 m b y 2 t e Y {x s u - x s d tan/3) 33 32 
+ eg tan/3) (l + e 3 tan/3) e// e// 



where 



e 3 = eg + y t 2 ey (9) 
a;f = (cos a, — sin a, z sin /3) (10) 
= (sin a, cos a, —i cos /3) (11) 

in the basis (S = H°,h°,A°). 

At large values of tan /3 the dominant supersymmetric contributions to rare decay B s — > 
H + H~ are mediated by neutral Higgs boson exchange that leads to [5j [25] 

BR(B - - 4 64 x 10" 6 M 2 ( ^ m ^^ \ " K^l) 3 T (u) 

bll{B s -+H ft )- 4.64 X 10 M Bs ^^2^27/4^3/21^1 J M 4 • i 12 ) 

Therefore, in this scenario, the magnitude of this observable is suppressed when \fiA t \ is 
small compared to Mg USY . 

As the gluino- quark- squark vertex is flavor diagonal for M ~ Msusy the dominant SUSY 
contributions to the b — > 57 rare decay come from the charged-Higgs boson and the chargino- 
stop loops. In particular the Wilson coefficients due to the charged Higgs contribution are 
proportional to the factor 



while the Wilson coefficient due to the chargino-stop loop has the form [28l [29 

\iA t tan (3 



7,S 



f( m h m h m l+) 



1 + 63 tan {3 

where / is the loop integral appearing at one loop. Eq 
contributions and the prescription used in Refs. [28| [29]. 



m b 



ht - dht 
m b 



1 + 63 tan (3 
2a, 



(14) 

includes the tan/3 resummed 

(15) 
(16) 



6ht = — t uM 3 tan j3 (cos 2 8- t C { 

37T 



mi L ,ml,Mi 



+ sin 2 etCoiml, m|, M 2 )) , 



(17) 



where, 5/i t is the correction to the charged-Higgs-top-strange vertex due to the gluino-stop 
loop and 9i is the stop mixing angle. 

The dominant supersymmetric contribution to the B u — > tv rare decay is due to the 
charged Higgs which interferes with the Standard Model contribution and we can define the 
ratio [30] 

-1 2 



R 



Btu 



B1Z(B U -> tv 
BK{B u ^tv) 



,MSSM 



SM 



1 - 



ni 



B 



m 



H± 



tan 2 /3 



1 + e tan (3 
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so as to quantify deviations from the Standard Model in this process. 

In addition, Ref. [21] has shown the importance of Kaon semi-leptonic decays in con- 
straining the charged Higgs contribution to the B u — > tv rare decay. In particular they 
consider the quantity 

K s (^ 2 )Kd(0+^0- 



V us (K l3 ) V ud [TTi2) 

where the subscript li refers to semileptonic decays with % final states and + — > + refers 
to beta decay. For the Standard Model, R123 = 1 while when a charged Higgs is included we 
have 



R 



1 



rn 



K 



rn 



H± 



1 



m d 



tan 2 /3 



1 + eo tan {3 



(20) 



The charged Higgs constribution in Eq. ([2"0l) and Eq. (|T8j) are the same and limits on R m 
can be a strong constraint on the scenario in which SUSY contributions to the B u — > tv 
dominate those of the Standard Model. Assuming that 5 is the largest allowed negative 
deviation of R123 from one and £ is the smallest allowed value of Rbtu, we see that for the 
charged Higgs to dominate over the SM contributions in Eq. (fT8|) the deviations must satisfy 
the constraint 



S > 



rn 



K+ 



m 



B u 



1 - 



m d 
m. 



j (1 + y/i) « 0.008(1 + y/i) 



(21) 



Hence, a two sigma experimental bound on 5 < 0.008(1 + would strongly disfavor 
scenarios in which the charged Higgs contribution to the B u — > rz/ decay is larger than that 
of the Standard Model. 
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2.1.2 M ~ Mqut 

When M ~ Mqut, corrections to the soft masses due to RG evolution are log enhanced. 
Therefore, if we neglect the yJYj term in Eq. (j2J), the left-handed down squark mass matrix 
is diagonalized by the matrix Ul which diagonalizes the up-quark mass matrix, rather than 
the down quark diagonalizing matrix Dl- Neglecting the corrections due to the bottom 
Yukawa over-estimates the splitting between the third and first two generations of down 
squark masses and is not valid when y^ ~ y t or when tan (3 is large. For \i and M 3 of the 
order of Msusy, with \xM 3 positive, one obtains e 3 ~ 0.01 and therefore the bottom Yukawa, 

m 6 tan/3 

Vb = -7— : 777 (22) 

v[l + e 3 tan p) 

becomes equal to yt for values of tan/5 > 100. The parametrization used in this article 
increases in accuracy as tan /3 takes smaller values, and also for larger values of fi, for which 
the above corrections to the bottom Yukawa coupling become significant, therefore reducing 
the value of y^. In this article, we shall assume that tan/3 < 60. In addition, as we shall 
discuss below, present experimental constraints lead to a preference for moderate or large 
values of /x at sizable values of tan/3 and small values of Ma- Therefore, we expect our 
parametrization to lead to a good approximation of the gluino induced effects in the scenarios 
discussed in this article. Furthermore using this approximation we were able to reproduce 
the the numerical B-physics limits obtained by Ref . [UJ , where the full renormalization group 
evolution of the mass parameters was performed. 

In the approximation, in which the left-handed down squarks are diagonalized by Ul, 
flavor violating vertices proportional to the CKM matrix in the gluino-down squark-down 
quark interaction vertex are induced, 

C g D V2g 3 ~g a (( V CK m) JI (d L ) J T a d L - (4)^4) , (23) 

and the soft SUSY breaking down-squark mass Lagrangian takes the form 

Cmass D (dlY(ml) I (d L ) J +(d* R ) I (ml) I (d R ) J + (24) 

where \x = /itan/3 — A&. Due to the gluino-quark-squark couplings being non-diagonal there 
are additional contributions to both the loop induced B s —>■ and b — > sj rare decays, 

but no large additional contributions to the B u — ► rv process. 

For M ~ Mqut, the effective flavor changing strange-bottom-neutral-Higgs coupling is [7] 

fY s sji _ ™dj(e 3 ~ e )(xg - x s d tan/?) t/3J , ta3J ( . 

{Xrl) "^(l + e tan/3)(l+ e3 tan/3) Ke//Ke// {2b) 

where we have assumed that the first two generations of left-handed down squark masses are 
mo, the uniform right-handed down squark soft mass parameters are m# and 

e w ^M 3 fiC (mlm 2 R ,\M 3 \ 2 ). (26) 
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In the limit of the left-handed sbottom mass being equal to that of the first two generations, 
Eq. (|25|) reduces to Eq. ([8]). In the M ~ M GUT scenario, the dominant SUSY contribution to 
B s — > // + //~ rare decay, at large tan /3, is found by substituting the form of X^ L in Eq. (l25l) 
into Eq. (fl2l) . The present experimental limit on B1Z(B S — > // + //~) disfavors very large 
positive contributions due to new physics effects. In high scale SUSY breaking, the new 
physics contributions to the B s — > fi + fi~ rare decay process are suppressed if the splitting 
in the left-handed down-type squarks soft mass parameters is such that e 3 — e is rendered 
small. As |eo| < |cq| this suppression may be significant whenever fiA t < 0, where the value 
of | //.At |, that allows such a cancellation, increases with the splitting of down squark masses 
and therefore with the messenger mass scale M. 

Futhermore, flavor violation in the gluino sector also induces relevant gluino contributions 
to the b — > S7 rare decay [MJ [35]. In Appendix [A] we find that, within the approximation of 
Eq. (123]) and Eq. ff24l) . the Wilson coefficients due to these gluino contributions are 

- g V2na 2 2 , M 3 e-"» / lAK (97 , 

7 ' 8 " Gf { ° QJ m b { ml (ml -m^) (ml -ml) {I7) 



m 2 bi (m 2 bi - ml) m 2 2 (m 2 2 - m 2 ) 

where cot 29 = (mg^ — m^J /(2\jl\mb), mQ 3 is the left-handed third generation down squark 
mass, is the i th sbottom mass, // = //tan/5 — At,, = arg(//), x g i = M 2 jm\. and the f^ g 
functions are defined in Eq. (|85|) q As expected this contribution to b — > 57 rare decay also 
disappears in the limit of uniform left-handed down squark soft mass parameters toq 3 = m$. 
For non-zero mass splittings, these contributions are important at large tan (3 and in the 
absence of CP violation are proportional to the sign of //M 3 . Therefore if //M 3 is positive the 
gluino contribution adds to that of the charged Higgs while when it is negative it subtracts 
from the charged Higgs contribution. 

In Fig. [1] we plot the relative gluino contribution to the decay branching ratio B1Z(b — > 
57), defined as 

e/toT-i/i w £5/v(0 > S7 J with gluinos O/vlO > -57/ without gluinos /oo\ 

8(BK(b -> S7)) = — fr r 2 , (28) 

LJ I'CyO ' 07 ) without gluinos 

for two different sets of SUSY parameters. The solid curve corresponds to the value of 
5(BlZ(b — > S7)) for superymmetric parameters // = M 3 = 300 GeV and X t = 0, while for 
the dashed curve we consider // = M 3 = —X t = 1 TeV. The splitting between the third and 
first two generations of squark masses is 20% and Ma ~ 200 GeV. We see that, in general, 
the gluinos lead to a moderate modification of BTZ(b — > 57). For instance, in the example 
shown in Fig. [fl the gluino effects lead to at most a 10-15% contribution to BlZ(b — > 57). In 
addition, since //M 3 is positive for these points, the gluino contribution to the B1Z(b — > 57) 
is also positive. 



A calculation of the gluino effects valid in the more general case has been recently performed in Ref. [36] 
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Figure 1: Variation of the gluino constribution to the b — > 57 rare decay branching ratio as 
a function of tan/3 for two different sets of SUSY parameters, assuming a 20% splitting in 
the squark masses and a value of Ma = 200 GeV. 



2.2 Direct dark matter detection through Higgs exchange 

The spin-independent elastic scattering cross-section for a neutralino scattering off a heavy 
nucleus is: 

4m? 



<?si 



7T 



(Zf p + (A-Z)f n 



where m r 



fp,n 



m^v is the mass of the nucleus, m x o is the neutralino mass, 



I J T q ^ 27 ^TG m „ 

\q=u,d,s q=c,b,t 



111 



p,n 



92m u 
4mwS{3 

92m d 
4m w C{3 



(g 2 N 12 - giN n ) 
(g 2 N 12 - giNn) 



1 



N\zS a c a I — - — 



h 



N u s a c a ( — 2 
V mi 



m 
1 



H 



+ JV14 Mr + 



m 



m. 



m H . 

0.118 ± 

0.062, /£ G « 0.84, = 0.014 ± 0.003, f$ d = 0.036 ± 0.008, /£ = 0.118 ± 0.062 and /£ G « 
0.83 [17]. In Eq. (I3T]) and Eq. (1321) A 7 ^- is the neutralino rotation matrix, a is the CP-even 
Higgs rotation angle and (m#) is the lighter (heavier) CP-even Higgs mass. Also in the 
above expression we define 



(29) 

(30) 

(31) 
(32) 



and the quark form factors are /£ 



0.020 ± 0.004, 



0.026 ± 0.005, /| 



1 + e tan p 
for the first two generations of quarks and 

m b 



1 + 63 tan (3 



(33) 



(34) 
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H.h 




Figure 2: Feynman diagram of the t-channel CP-even Higgs contribution to the spin- 
independent cross-section. 

for the bottom quark. In Eq. (1321) . we are ignoring the contribution from s-channel squark 
exchange, which becomes subdominant for heavy squark masses. In the limit of large tan (3 
ad is tan/3 enhanced compared to a u . Moreover, for large tan/?, \i 3> Mi, M 2 ~ M\ and 
120 GeV < M A < 600 GeV (M A < 120 GeV), one obtains N u > N 12 , m H ~ M A 
(rrih — Ma) and s a ~ — l/tan/3 (c a ~ l/tan/3). Hence we find that the dominant 
contribution is 

f „ _ m ( K±K i 2 /rg V^iVniV^tan/? 

Jp ' n ~ p ' n Vl + e tan /? + 27 1 + e 3 tan (3 ) Am w M\ 1 } 

» f °' 14 | °-° 6 \ m 9i92NnN 13 tanP 
V 1 + e tan /? 1 + e 3 tan /3 / p 4m^M] 

where we have neglected, in the first line, the splitting between the first two generations of 
squarks and, in the second line, the differences between the proton and the neutron mass and 
we used the fact that the neutron and proton fx factors are relatively similar. Assuming that 
the mass of the neutralino is much larger than that of the nucleus we have m r ~ ~ Am p 
and 



AA 2 m 2 p 

°si glolNlxN^ml tan 2 (3 ( 0.14 0.06 



osi « P -A 2 f v (37) 



A 4 ' A-nm 2 w M\ \ 1 + e tan (3 1 + e 3 tan (3 ) ' K ' 

where asi/A 4 is the neutralino nucleon spin-independent cross-section. From Eq. (|38|) the 
spin-independent cross-section scales as tan 2 f3/M\ and therefore direct dark matter detec- 
tion experiments like CDMS can put strong constraints on regions of small M A and large 
tan f3. 
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3 Numerical limits and constraints 

3.1 Experimental constraints on B-physics observables. 

Due to the extra supersymmetric contributions to these rare decays we find that experimental 
data put strong constraints on these models. In particular, the world experimental average 
of the branching ratio of the b —>■ s'j rare decay is [38] 

BK(b -> S1 ) exp = (3.55 ± OMt°o°l ± °- 03 ) x 10 ~ 4 > ( 39 ) 

which agrees well with the Standard Model prediction [39] 

BK(b -> s 7 ) 5M = (3.15 ± 0.23) x 10~ 4 . (40) 

Using the experimental and SM ranges for the b — > sj rare decay we find the 2a allowed 
range is 

BV(b -> s >y) MSSM 
0.89<iW,= skw^ <41) 

For the B u — > tv rare decay the Belle experimental collaboration measures a branching 
ratio of HQI 



BK{B U - tv) Bc11c = (1.79l°;l(stat)t°;^(syst)) x 10~ 4 , (42) 

while the Babar collaboration finds the preliminary value [UJ 

BTZ{B U -> rz/) Babar = (1.20 ± 0.54) x 10" 4 . (43) 

The average of these two experiments is then [42] 

Bft(.B u -> rz/) Exp = (1.41 ± 0.43) x KT 4 . (44) 

Using / fl = 189±27 MeV from LQCD [42] and the average value of \V ub \ = (3.98±0.45) x 10~ 4 
from HFAG [43J, the Standard Model prediction is 

BK{B U -> tv) = (1.09 ± 0.40) x 10~ 4 (45) 

Assuming at most a 2cr deviation from new physics, we find the allowed range 

B1Z( B —, t1j \ m ssm 

om * - Zk £ 25L (46) 

For the R123 constraint in Eq. (I20l) . Ref. [3TJ finds that 

0.990 < R m < 1.018 (47) 
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when they use the value of fx/ fn = 1-189 ± 0.007 from Ref. [32]. However if we use the 
average value of fx/ fn = 1.19±0.015 from Ref. [33] 

0.96 < ife < 1-05. (48) 

In this paper we will use this more conservative limit on the -R/23 rather than the one in 
Eq. (1471) . From Eq. ( J2TT) it is clear that the more restrictive bound in Eq. ( 1471) strongly 
disfavors the regions of low Ma and large tan f3 that are allowed by the B u — > tv constraint 
in Eq. (I4"6"]) . where the charged Higgs contribution to the B u — > tv process dominates that 
of the Standard Model. 

The B s — > rare decay has yet to be experimentally observed. The present experi- 

mental exclusion limit at 95% C.L. from CDF [44] is 

BU{B S -> < 5.8 x 10~ 8 , (49) 

which can put strong restrictions on possible flavor changing neutral currents in the MSSM 
at large tan/3. Additionally the projected exclusion limit, at 95% C.L., on this process for 
4 fb _1 at the Tevatron is [45] 

BU{B S -> < 2.8 x 10" 8 . (50) 

For the LHC, the projected ATLAS bound at 10 fb" 1 is [46] 

BU{B S -> /i + /i~) < 5.5 x 10" 9 . (51) 

In addition, LHCb has the potential to claim a 3a (5a) evidence (discovery) of a standard 
model signature with as little as ~ 2fb~ 1 (6fb~ 1 ) of data [47]. 

3.2 Direct dark matter detection constraints. 

As discussed in Section 12.21 the spin-independent scattering cross-section for a neutralino 
off a heavy nucleus scales as tan 2 /3/M^ and therefore puts strong constraints on the SUSY 
parameter space. At present the CDMS [20] and XENON [48J collaborations have put a limit 
on the spin-independent neutralino-nucleon cross-section that is of the order of 10~ 7 pb [21 J . 
We will use the CDMS limits throughout this paper because we only consider neutralino 
masses greater than 100 GeV and the current limits from the CDMS experiment are slightly 
stronger than those from XENON [21] for these range of masses. By the end of 2009, the 
sensitivity of the XENON100 experiment will improve by an order of magnitude [22], while 
the sensitivity of the SuperCDMS experiment will improve by factor of 5 [23]. Therefore in 
the near future, these direct dark matter detection experiments will be able to probe regions 
of SUSY parameter space that will also be probed by the Tevatron in non-standard Higgs 
searches. 
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3.3 Parametric scenarios 



Due to the dependence on the messenger scale M we consider different scenarios to illustrate 
the interplay between B-physics, Higgs physics and dark matter searches within the frame- 
work of the MSSM with large tan /3. In this section we will assume the gaugino unification 
relation |M 2 | ~ 2|Mi| and that all the gauginos have equal phases. We shall consider four 
parameters: the CP-odd Higgs mass M4, the ratio of the two Higgs vacuum expectation val- 
ues tan /?, the Higgsino mass parameter \i and X t = A t —jj,/ tan /3, where A t is the stop-Higgs 
trilinear coupling. In order to consider the direct dark matter detection constraints, we will 
also assume that the total relic density agrees with that found by WMAP, independently 
of the squark, neutralino and Higgs spectrum. This may require an appropriate slepton 
spectrum or a departure from the standard thermal dark matter predictions [H]. In addi- 
tion we considered the non-standard Higgs boson search limits in Ref. [241 [50| |5T] and used 
the CPsuperH [52] program to project these constraints onto the Ma — tan (3 plane. Even 
though the Tevatron will be collecting about 8 fb _1 of data, the projected 4 fb _1 CDF limit 
provides a conservative estimate of the final reach of the Tevatron in the H/ A — > rr channel 
because a realistic treatment of the detector and efficiencies may lead to somewhat weaker 
constraints than those shown in Ref. [50J. The future LHC constraints on the H/A — ► rr 
channel, from Ref. [51], correspond to the projected limits at ATLAS for 30 fb _1 of data. 

The good agreement between the Standard Model prediction and the experimental mea- 
surement of the branching ratio of b — > 57 implies that either there is some cancellation 
between the dominant new physics leading order Wilson coefficients in Eq. (1131) and Eq. (TT41) 
or each of them are individually small. If A t is sizable and the sign of \xA t is negative then 
a cancellation between the chargino-stop and charged Higgs Wilson coefficients is possible 
for large enough values of tan (3. For negligible values of A t the chargino-stop contribution 
is suppressed, which requires the charged Higgs amplitude to be small. This suppression of 
the charged Higgs Wilson coefficient can be achieved by making the bottom-top-charged- 
Higgs vertex small through a cancellation between the tree-level coupling and the one loop 
tan/? enhanced correction in Eq. ffTBl . In addition to the leading order contribution we have 
also included the next-to-leading order contributions due to the charged Higgs as discussed 
in Ref. [37J. Additionally there is also a possible gluino contribution that can be signifi- 
cant if the squark masses are splitted. The gluino contribution, for the scenario in which 
M ~ Mqut, is relevant for large //, M3 and tan/3 and depending on the sign of ^3 this 
contribution interferes constructively or destructively with the other two contributions. 

Similarly to the b — > 57 rare decay, the B s — > fi + fi~ rare decay is sensitive to the scale 
M at which supersymmetry breaking is communicated to the visible sector. Depending on 
the scale M the structure of the X^ L couplings in Eq. (jSJ) and Eq. (J2"5l) is different. If 
M ~ Msusy, this coupling is suppressed if ey or equivalently A t is small compared to m or 
M3. However if M is large compared to Msusy then splittings in the squark masses, due to 
RG evolution, can induce a suppression of this coupling, due to a cancellation between the 
stop-chargino and sbottom gluino loops. 

The B u — ► tv bound in Eq. (|46p imposes a complementary constraint to that of b — > 37. 
The lower bound of Rb T u implies that there cannot be complete destructive interference 
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Figure 3: Plot of X t versus /x, where M4 = 110 GeV, tan/3 = 40, sign(X t ) negative, 
M3 = 800 GeV. The green (grey) hatched region is allowed by the B u — > rz/, 6 — > 37 and 
-S« - ► /x + /x~ constraints for M ~ Ms;ysy while the yellow (light grey) region is allowed by 
the same constraints for the M ~ Mqut scenario. The region below the horizontal red (dark 
grey) lines has been probed by the CDMS direct dark matter detection experiment assuming 
that the LSP is mainly bino and \M\\ = 2\M 2 \. The solid (dashed) line corresponds to a 
Wino mass parameter M 2 = 200 (500) GeV. 

between the SUSY and Standard Model contributions. Therefore there are two disconnected 
allowed regions in the Ma — tan /3 plane: one where the charged Higgs induced amplitude 
dominates the Standard Model contribution and the other where the opposite happens. 

In these kind of scenarios, Eq. (|38|) suggests that the spin- independent dark matter scat- 
tering cross-section is quite sensitive to the amount of the Higgsino component in the lightest 
neutralino. At low values of /x there is a large Higgsino component to the lightest neutralino 
and hence a larger scattering cross-section through t-channel CP-even Higgs bosons. For 
large values of /x, instead, the Higgsino component is much smaller, and so the coupling of 
the neutralino to the Higgs is suppressed leading to a smaller cross-section. In addition, the 
sensitivity of direct dark matter detection experiments, like CDMS, to the spin independent 
cross-section decrease with increasing LSP mass for m Xl ^ 50 GeV. Observe, however that 
the CDMS constraint assumes the neutralino to be the dark matter candidate while in low 
SUSY breaking scenarios the LSP is naturally the gravitino, and therefore these constraints 
should not apply. 

Throughout this section we set the uniform left-handed soft squark mass parameter for 
the first two generations to m = 1.5 TeV, the third generation soft squark mass parameters 
are mu a = mq z = 1.2 TeV and the uniform right-handed down squark soft mass parameter 
is mji = 1.5 TeV. This form of the left handed down squark soft masses has been chosen so 
as to mirror a 20% splitting in the soft masses due to their renormalization group evolution 
in M ~ Mqut scenario, as is naturally the case whenever the gaugino masses are of the same 
order as the scalar masses at the scale M. We have chosen a value of the third generation 
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Figure 4: Plot of X t versus /i, where Ma = 200 GeV, tan/5 = 55, sign(X t ) negative, 
M3 = 800 GeV. The green (grey) hatched region is allowed by the B u — > ru, b — > 57 and 
B u ^ + ^ constraints for M ~ Msusy while the yellow (light grey) region is allowed by 
the same constraints for the M ~ Mqut scenario. The region below the horizontal red (dark 
grey) lines has been probed by the CDMS direct dark matter detection experiment assuming 
that the LSP is mainly bino and \M\\ = 2\M2\. The solid (dashed) line corresponds to a 
Wino mass parameter M 2 = 200 (500) GeV. 

squark masses slightly larger than 1 TeV, in order to satisfy the lightest CP-even Higgs mass 
constraints for all the scenarios under study. For the M ~ Msusy scenario we shall assume 
a degenerate squark spectrum with soft masses of 1.2 TeV. In this way, we can compare 
the results with those in the case M = Mqut for which the third generation masses, most 
relevant in the calculation of the B-physics observables, have also values of 1.2 TeV. 

In Fig. [3] and Fig. H] we present the effects of the B physics constraints and the CDMS 
direct dark matter detection experiment limit on the X t — \x plane for (Ma, tan (3) = 
(110 GeV, 40), and (Ma, tan P) = (200 GeV, 55) respectively. These two sets of values 
correspond to regions of parameter space which are close to being probed at the Tevatron 
in inclusive A/H — > rr searches [23]. The green (grey) hatched region is the one allowed by 
the b — > sj, B s — > fi + n~ and B u — ■> rv constraints for M ~ Msusy, while the yellow (light 
grey) region is allowed by the same constraints for M ~ Mqut scenario. The B u — > rv 
constraint does not depend on the parameter X t and therefore the constraint in Eq. fj46l) 
selects a horizontal band in the X t — \x plane. The regions below the solid and dashed red 
(dark grey) lines is excluded by CDMS, for M 2 = 200 GeV and M 2 = 500 GeV, respectively. 

In Fig. [31 for M4 = 110 GeV and tan/3 = 40, The extra gluino contributions to the 
b — > S7 and, most relevantly, the B s — > fi + n~ rare decay rates, leads to a modification of 
the preferred values of X t . While the B-physics constraints lead to a preference for small 
values of X t in the M ~ Msusy scenario, moderate values of X t ~ —500 GeV are preferred 
in the M ~ Mqut case. Assuming that the LSP is the neutralino, which is natural in the 
M ~ Mqut scenario, the recent CDMS limits [20j are quite strong and exclude regions below 
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Figure 5: Plot of M A versus tan 0, where X t = -400 GeV, y, = 800 GeV and M 3 = 800 GeV. 
The green (grey) hatched region is allowed by the B u — ► tu, b — *■ S7 and I? u — > y + y~ 
constraints for M ~ Ms;ysy while the yellow (light grey) region is allowed by the same 
constraints for the M ~ Mqut scenario. The region above the red (dark grey) lines has 
been probed by the CDMS direct dark matter detection experiment assuming that the LSP 
is mainly bino and \M\\ = 2\M2\. The blue-green region is excluded by Non-standard Higgs 
boson searches in the inclusive rr channel at 1.8fb _1 . The region above the black solid 
(dashed) lines will be probed in the H/A — > rr channel at the Tevatron (LHC) with a 
luminosity of 4 fb _1 (30 fb" 1 ). 

~ 900 GeV. Indeed, for these values of Ma and tan/3 we observe that large to moderate 
values of \i are preferred for both SUSY breaking scenarios. 

Fig. H] shows the situation for Ma = 200 GeV and tan/3 = 55. Similar to Fig. [3] small 
values of X t are preferred in the M ~ Msusy scenario while moderate values of X t ~ 
—400 GeV are preferred in the M ~ Mqut scenario. For these values of Ma and tan/?, 
the CDMS experimental bound is less stringent than in Fig. [31 restricting values below 

~ 500 GeV in the M ~ M GUT scenario. Similar to Fig. [3], the B s — > fi + n~ constraint is 
the main discriminant between the M = Mqut and M = Msusy scenarios. 

From Figs. [3] and Fig. H] we can observe some generic features. Independent of the 
SUSY breaking scale, for these regions of parameter space, that can be probed at the Teva- 
tron, one obtains that X t < 0.5Msusy- These low values of the stop-mixing parameter 
X t imply an upper bound on the lightest CP-even Higgs boson mass, m h < 120 GeV and 
therefore could be within the reach of the Tevatron collider. As we had previously em- 
phasized the regions close to (Ma, tan P) = (110 GeV, 40) and (200 GeV, 55) are yet to be 
probed at the Tevatron in H/A — ► rr searches at 1.8 fb _1 . Using Figs. [3] and H] we see 
that the region of parameter space around (Ma, tan (3, X t , y) ~ (110 GeV, 40,0, 1 TeV) and 
(M A , tan 0, X t , y) ~ (200 GeV, 55,0, 1 TeV) satisfies all the constraints in the M ~ M SU sy 
scenario. In addition, from Fig. HJ we also find that the region of parameter space close 
to (M A Mn(3,X t ,n) ~ (200 GeV, 55, -400 GeV, 800 GeV) satisfies all constraints for M ~ 
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Figure 6: Plot of M A versus tan(/3), for X t = 0, n = 1000 GeV and M 3 = 800 GeV. The 
green (grey) hatched region is allowed by the B u — > tv, b —>■ sj and B u — > fi + n~ constraints 
for M ~ Msusy while the yellow (light grey) region is allowed by the same constraints for 
the M ~ Mqut scenario. The region above the red (dark grey) lines has been probed by 
the CDMS direct dark matter detection experiment assuming that the LSP is mainly bino 
and | Mi | = 2|M2|. The blue-green region is excluded by Non-standard Higgs boson searches 
in the inclusive rr channel at 1.8fb _1 .The region above the black solid (dashed) lines will 
be probed in the H/A — > rr channel at the Tevatron (LHC) with a luminosity of 4 fb -1 
(30 fb- 1 ). 

Mqut- As the constraints from B-physics, Higgs physics and direct dark matter searches get 
stronger this kind of analysis could help us identify regions of parameter space that would 
still be compatible with all experimental limits. 

In Fig. we consider the /i = 800 GeV and X t = —400 GeV parametric scenario. 
The region above the solid and dashed red (dark grey) lines has been probed by CDMS in 
direct dark matter detection experiments, for a Wino mass parameter M 2 = 200 GeV and 
M 2 = 500 GeV respectively. The blue-green (medium grey) region is excluded by CDF and 
DO in non-standard Higgs boson searches in the rr channel at 1.8 fb -1 . The green (grey) 
hatched region is allowed by the experimental constraints on the b —>■ 37, B s — > and 
B u — > tv rare B decays for M ~ Msusy while the yellow (light grey) region corresponds to 
the same constraints for the M ~ Mqut scenario. The region above the black solid (dashed) 
lines will be probed in the H/A — > tt channel at the Tevatron (LHC) with a luminosity of 
4 fb -1 (30 fb -1 ). In this parametric scenario, with \iX t < and /iM 3 > and the different 
b — ► 57 contributions tend to cancel against each other. The chargino-stop approximately 
cancels the charged Higgs contribution in the M ~ Msusy scenario, while in the M ~ Mqut 
scenario the chargino contribution tends to cancel both the charged Higgs and the gluino 
contributions. There are two regions allowed by B u — > tv constraint. The allowed region 
at low Ma and large tan/5 is where the supersymmetric contribution dominates, while in 
the large Ma and low to moderate tan j3 region the Standard Model contribution dominates. 
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The B s — > yU+yU - constraint is quite strong for the M ~ Mj^sy scenario due to a non-zero 
X t , excluding most of the region where the charged Higgs contribution to B u — > rz/ becomes 
dominant, while for the M ~ Mqut a cancellation is induced in the flavor violating effects 
due to the splitting of the left-handed down squark soft mass parameters. 

In Fig. [5] we see that somewhat smaller values of tan (3 are allowed in the M ~ Msusy 
case than in the M ~ Mqut case. This effect may be explained by the gluino contributions 
to b — > 57 : For positive values of fiM 3 a larger chargino-stop contribution is required to get 
agreement with the experimental values, which may be obtained for larger values of tan/3. 
The R123 constraint in Eq. (j4"5j) becomes too weak to give any significant constraint in Fig. [5] 
and in the other scenarios we consider in this paper. Let us stress again, if we had considered 
the more restrictive bound on i?; 2 3 in Eq- dSD the region of low Ma and large tan f3 allowed 
by the other flavor constraints would be strongly disfavored. 

As Fig. [3] suggested the region of small Ma and tan/3 ~ 35-45 is allowed by all the B 
physics constraints and has yet to be probed by the Tevatron in inclusive A — > tt searches. 
The allowed region is larger for M ~ Mqut than for M ~ Msusy- For M ~ Mqut, this 
region is however also constrained by direct dark matter detection experiments like CDMS 
as can also be seen in Fig. [51 In addition, there is also a region around (M^,tan/3) = 
(200 GeV, 55) that is allowed in the M ~ Mqut scenario, which has yet to be probed in 
direct dark matter detection experiments and non-standard Higgs searches. However the 
XENON100 and SuperCDMS experiments should be able to probe this region of parameter 
space in the near future due to their improved sensistivity. 

In Fig. [6] we consider a scenario where /x = 1 TeV and X t = 0. The chargino-stop 
contribution to the b — > 57 rate is small in this scenario, because At ~ 0, while the charged 
Higgs contribution tends to be suppressed because of a cancellation between the 1-loop 
and 2- loop contributions in Eq. (1131) . The B s — > fi + fi~ constraint for the M ~ Mqut 
scenario is strong because there is no cancellation that occurs in Eq. (1251) . while it is weak 
in the M ~ Msusy scenario because At ~ 0. In addition there are two regions around 
(M A , tan P) = (175 GeV, 55) and (M A ,tan/3) = (115 GeV, 40) that are allowed by all these 
constraints in the M ~ Msusy scenario but disallowed in the M ~ Mqut scenario. As in 
the previous case, a significant region of parameters consistent with all experimental appears 
at low values of Ma and large values of tan f3 but in this case is compatible with M ~ Msusy 
and hence no CDMS restrictions apply. In addition both Fig. [5] and Fig. [6] also show that the 
Tevatron collider will be able to probe all these allowed regions in non-standard Higgs boson 
searches with a luminosity of 4 fb _1 . Futhermore, in both these scenarios, the LHC will be 
able to probe more of the allowed regions of large Ma and low tan/3, in the H/A — > tt 
channel with 30fb _1 of luminosity. 

4 Conclusion 

In this article we have studied the effect of varying the messenger scale on B physics observ- 
ables within Minimal Flavor Violating supersymmetric models. In particular we found that 
the b — > 57 and B s — > rare decays are sensitive to the scale M at which supersymmetry 
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breaking is communicated to the visible sector. Considering the effects of the RG evolution 
which amounts to an alignment of the left-handed squark masses with the up Yukawa cou- 
plings, with uniform right-handed down squark soft masses and also uniform left-handed 
down squark masses of the first two generations, we have derived an analytic expression for 
the gluino contribution to the b — > rare decay. We find that the gluino contribution is 
proportional to the splitting between the third generation left-handed down squark mass 
and that of the first two generations. The relative sign of the gluino contribution to that 
of charged Higgs depends on the sign of /iM 3 . Hence in the case of the messenger scale 
M ~ Mqut, when the splitting in the left-handed squark masses is non-zero, this contri- 
bution can be significant. In addition we also show the dependence of the dominant SUSY 
penguin contributions to the B s — > fi + fi~ rare decay on the scale M. 

We have also studied the interplay between the B-physics constraints, dark matter direct 
detection experiments and non-standard Higgs boson searches at the Tevatron. For large 
soft squark masses, the spin-independent neutralino nucleon cross-section is proportional to 
tan 2 f3/M\ and hence direct detection experiments put strong constraints on regions of low 
M A (M A < 200 GeV) and large tan/3. In particular, we have projected the CDMS direct 
dark matter detection experimental constraint on the Ma — tan (3 plane, for different values 
of n and M 2 . Independently of the messenger scale M, the B physics, Higgs physics and 
Dark Matter experimental constraints suggest that low values of X t and large to moderate 
values of /i are preferred. Such low values of X t generally suggests an approximate upper 
bound on the lightest Higgs mass m h < 120 GeV, which may be within the reach of the 
Tevatron collider. 

In addition, we have presented parametric scenarios that satisfy all the B physics experi- 
mental constraints considered in this article, within the scenarios of low scale (M ~ Msusy) 
and high scale (M ~ M GUT ) supersymmetry breaking, and can be probed by the Tevatron 
collider and the LHC in the near future. In particular for M ~ M GUT we find a region 
around (Ma — 200 GeV and tan/3 ~ 55, and moderate values of \X t \ and fi, which is within 
the 4fb _1 reach of the Tevatron collider in non-standard Higgs searches. For M ~ Msusy, 
instead, smaller values of \X t \ and moderate or large values of /i are preferred, in order to 
obtain acceptable values of B1Z(B S — > and B1Z(b — > s^y). Moreover, we showed that 

for X t ~ 0, there are large regions of parameter space for low M A and large tan (3 that remain 
to be probed by non-standard Higgs searches at the Tevatron collider. Apart from a region 
at similar values of Ma and tan f3 as the ones arising in the M ~ Mqut scenario, we found 
an additional region, for smaller values of M A ~ 115 GeV and tan/5 ~ 40. This region 
only appears in scenarios with low energy supersymmetry breaking, since it is constrained 
by direct dark matter searches in the M ~ M GUT scenario. 

Our analysis suggests that, in Minimal Flavor Violating MSSM, future Higgs searches at 
the Tevatron and direct dark matter detection experiments at CDMS and XENON will reveal 
useful information about SUSY parameters and the scale of supersymmetry breaking. For 
instance, the detection of a light non-standard Higgs boson at the Tevatron and dark matter 
at the XENON and CDMS experiments, by the end of 2009, would suggest a M ~ M GUT 
scenario, from which we can infer moderate values of X t , large to moderate values of ji 
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and a Standard Model Higgs boson mass < 120 GeV. On the other hand, detection of a 
light non-standard Higgs boson and non-detection of dark matter may suggest a lower SUSY 
breaking messenger scale, small values of X t , large values of /i and a Standard Model Higgs 
boson mass close to that of the LEP experimental limit. 
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A Gluino contribution to b —> 57 

In the initial gauge basis the gluon-quark-squark interaction Lagrangian has the form 

C 9 D V2g 3 ~g a (K) W L - (4)^4) . (52) 
Rotating the quarks into the mass basis by the matrices 

di->Dg"d£ 4-^4 (54) 
and the down squarks by the matrices 

dl^Ul J d J L 4^^ J 4 (55) 

so as to diagonalize the down squark soft masses. Hence the gluon-quark-squark interaction 
Lagrangian becomes 

C g D V2gfg a {{uId L ) J 1 {d*^ T a d[ - (4)^4) (56) 

= V2g 3 g a (jy C KM) JI (dr L ) J T a d I L - (4)^4) • (57) 

However this rotation induces off-diagonal terms to the left-right and right-left blocks of the 
down-squarks, so that, 

C mass D (dl) I (miy(d L ) I + (d* R ) I (ml) I (d R ) I + (58) 

where we have neglected terms proportional to the down and strange Yukawa's and fx = 
/itan (3 — A b . Hence all three of the left-handed states mix with rn 3 dR . As V^ 3 KM xs 1 we 
explicitly diagonalize the {d\, d R ) sector so that 

4 -> ceh-e'^seh (59) 
4 e^seh + cek (60) 
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where 
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To leading order in the CKM matrix elements we can further diagonalize the left-handed 



7 = 1,2 states so that 
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Hence 
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where we have neglected higher orders in the CKM matrix. Therefore the vertex factors are 
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Now rewriting the following identity in Eq. (Jb! 

c s e (m^ - mjy + |//|m 6 (cg - s 2 e ) 



we have 
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Hence using the amplitudes defined in Ref. pj we find the Wilson coefficients 
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where 



= 9^ + W L T ? ]ogx (85) 
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